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We calculated spectra of the first six members of the Mg-like isoelectronic sequence using different 
approximations. The most accurate results were obtained with the configuration interaction -|- all 
order method (CI-I-AO), which provided relative accuracy for transition energies on the level of 
0.1%, or better. On this level of accuracy the Breit and QED corrections become important for 
the systems with nuclear charge Z > 20. The retardation part of the Breit interaction and QED 
corrections for the partial waves with I A 0 are still too small to be important. 

PACS numbers: 31.15.ac, 31.30.jf. 


I. INTRODUCTION 

High accuracy atomic calculations are necessary not 
only for atomic physics itself, but also for different appli¬ 
cations from atomic frequency standards to tests of fun¬ 
damental symmetries, search for the variation of the fun¬ 
damental constants, and astrophysics. Atomic Mg and 
ions of the Mg isoelectronic sequence are convenient sys¬ 
tems to test theoretical methods for atomic calculations 
mu- With ten core and two valence electrons they 
require accurate treatment of the core, core-valence, and 
valence correlations. Relativistic and quantum electrody¬ 
namic (QED) corrections are very small for magnesium, 
but grow along the isoelectronic series. Therefore in the 
high accuracy calculations of the heavier ions we need to 
account for both electronic correlations and relativistic 
and QED corrections. 

At present there are several methods to calculate spec¬ 
tra and other properties of the many-electron atoms. The 
many-body perturbation theory (MBPT) is quite effec¬ 
tive for monovalent atoms [mis]. However, the order- 
by-order approach leads to rather slow convergence and 
becomes impractical above the third, or the fourth or¬ 
der. Because of that different variants of the all-order 
(AO) methods are currently used instead [T3]. The cou¬ 
pled cluster (CC) method is the most popular all-order 
method, which is used by several groups [TKH2H] . 

The MBPT-based approaches, including the CC 
method are most effective in treating core and core¬ 
valence correlations, but are less suitable for treating va¬ 
lence correlations, where there is no well defined small 
parameter. Here the non-perturbative approaches, such 
as configuration interaction (Cl) [21], multi-configuration 
Hartree-Fock (MCHF), or multi-configuration Dirac- 
Hartree-Fock (MCDHF) |2l II] are more effective. They 
allow accurate treatment of correlations between several 
valence electrons, but start to fail when the number of 
correlated electrons exceeds four, or five. 

There are also mixed approaches where core and core¬ 
valence correlations are treated perturbatively and va¬ 
lence correlations are treated within Cl method. The 
simplest method of this type includes Cl and the second 
order MBPT (CI-I-MBPT) [T1|S1[25|. The more advanced 


variant includes Cl and linerized CC (CH-AO) [H ITT] . 
In both cases we use either MBPT, or CC method to 
form an effective Hamiltonian for the valence electrons 
and then we use Cl method to find valence energies and 
many-electron valence wave functions. 

Relativistic effects for many-electron atoms are usually 
included within Dirac-Coulomb, or Dirac-Coulomb-Breit 
no-pair approximation. Sometimes, one-electron QED 
corrections (Lamb shifts) are also included using effec¬ 
tive, or model potentials [261130] . 

In this paper we study relative size of different correc¬ 
tions to the transition energies in Mg and Mg-like ions 
up to Cl VI. We use CH-MBPT and CH-AO methods 
for Dirac-Coulomb and Dirac-Coulomb-Breit Hamiltoni¬ 
ans mm- Lamb shift corrections are included only for 
the s-electrons. We find that the relative accuracy of 
the CI-I-AO method for the transition frequencies is on 
the order of 0.1% and remains almost constant along the 
sequence. At the same time, the dominant theoretical 
errors are different in the beginning and in the end of the 
sequence. While the role of the higher order correlation 
corrections is decreasing, the role of the relativistic and 
QED corrections grows. 


II. THEORY AND METHOD 

In this paper we use variants of the CI-I-MBPT and 
CI-I-AO methods, which are based on the Brillouin- 
Wigner perturbation theory. The resultant effective 
Hamiltonian for the valence electrons is symmetric, but 
energy-dependent mm- Savukov and Johnson sug¬ 
gested an alternative approach based on the Rayleigh- 
Schrodinger perturbation theory. In this case the effec¬ 
tive Hamiltonian is non-symmetric and does not depend 
on the energy. 

We do not use variant of Savukov and Johnson because 
of the well-known problem of intruder states. In our ap¬ 
proach we do not have intruder states, but we need, in 
principle, to account for the energy-dependence of the ef¬ 
fective Hamiltonian. Fortunately, if one is interested in 
the low-lying atomic states and forms effective Hamilto¬ 
nian according to the recipe from Ref. [T] , the respective 
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corrections are very small. Thus, in the first approxima¬ 
tion one can calculate MBPT diagrams at fixed energies. 
These energies are chosen to be the Hartree-Fock energies 
of the first valence orbitals with given symmetry (i.e. for 
Mg-like ions the energy £ 3 ^ is used for all valence s elec¬ 
trons, etc.) [ 1 ]. In the next approximation one can also 
calculate first derivatives of all diagrams in respect to 
the energy. The effective Hamiltonian for a given energy 
is then formed by extrapolating diagrams to this energy 
[nimisi]. Here we use this method to study respective 
corrections to the theoretical spectrum. We conclude, 
that at the present level of accuracy these corrections do 
not improve agreement with the experiment and can be 
neglected. 

Our Cl and CI-I-MBPT calculations are done with 
the package described in Ref. [32]. CI-I-AO calculations 
are done with the extended variant of the same pack¬ 
age. These calculations include several steps. At first 
we solve linearized coupled cluster equations for the core 
and core-valence cluster amplitudes in the single-double 
(SD) approximation. Then we form one-electron and 
two-electron valence cluster amplitudes. On the third 
step we use these amplitudes to form effective valence 
Hamiltonian (see Ref. El for details). In our opinion the 
CI-I-AO method described above treats valence-valence, 
valence-core, and core-core correlations in a most effec¬ 
tive manner. 


Calculation details 


eluded in the Cl space, in the pure Cl calculation we ne¬ 
glect all core-valence and core-core correlations. In the 
CI-I-MBPT and CI-I-AO calculations these correlations 
are included in the effective Hamiltonian of the valence 
electrons. 


QED corrections 

We use semiempirical approach to account for QED 
corrections where we include only the Lamb shift for the 
s-electrons [It is known that Lamb shift for other partial 
waves is at least one order of magnitude smaller]. Follow¬ 
ing |30j we can parametrize Lamb shift for the hydrogen¬ 
like ion as (we use atomic units h = e = rrie = 1 ): 

Aeqed = - irFn{aZ ), (I) 

Trn'^ 

where n is the principle quantum number, a is the fine 
structure constant, and Z is the nuclear charge. This 
expression can be generalized for the non-diagonal matrix 
elements: 

a^Z'^ 

<^n,n' = 7 A3/2 ' (^) 

Function F weakly depends on aZ and on indexes n,n'. 
This function is tabulated in Ref. m- Noting that 
^(„^q 3/2 = |«'n(0)4'n'(0)| we get: 


We calculate spectrum of each ion of the isoelectronic 
sequence within several approximations, namely pure va¬ 
lence Cl, CI-I-MBPT, and CI-fAO methods. We use 
Dirac-Coulomb Hamiltonian in the no-pair approxima¬ 
tion. We also do calculations with Breit and QED correc¬ 
tions. Finally, for the CI-I-MBPT and CI-I-AO we repeat 
computations with account for the energy dependence of 
the effective Hamiltonian. Previously this dependence 
has been studied for the CH-MBPT [TJ |3T] , but not for 
the CI-I-AO method. At the same time, the consistent 
treatment of the high orders requires accurate account of 
the energy dependence of the effective Hamiltonian [13] . 

For each ion we form the finite basis set, which in¬ 
cludes Dirac-Fock orbitals for the core and valence states 
and Sturmian orbitals for virtual states. Sturmian or¬ 
bitals effectively account for both discrete and continu¬ 
ous spectrum [35]. After adding Sturmian orbitals we 
diagonalize Dirac-Fock operator of the core on the whole 
basis set (i.e. we use the V^~‘^ potential, where = 12 
is the total number of electrons). 

This way we make 30spdfgh basis set of 300 orbitals 
that includes partial waves with Z = 0,..., 5. We use the 
whole basis set for the MBPT and CC parts of the calcu¬ 
lation. The configuration space for the Cl part includes 
single (S) and double (D) excitations to all orbitals up 
to 21spdfg. This corresponds to the full two-electron Cl 
on the basis set 21spdfg. As no core excitations are in- 


a'^Z 


^n,n' — 


^'n(0)«'„-(0)|F„,„,(aZ). 


(3) 


This expression can be used not only for hydrogen-like 
ions, but also for the valence electrons of many-electron 
atoms. Following Ref. [3S] we can express electron den¬ 
sity at the origin in terms of the binding energy 


l'I'.( 0 )'I'.'( 0 )| 


jz, + lyz 

+ 1 ) 


(4) 


In these expressions Zi is the charge of the ion (for neutral 
atom Zi = 0 ) and v is the effective quantum number. 
Now we can write QED corrections as: 

_ (4£ j /£ j /') ^ „ 3^7-2 771/' „ ry \ (^\ 

^v,v’ — I 1 ', OL Z F{aZ). (5) 

7r(Zi -I- I) 


Binding energy Si, for the valence electrons in the many- 
electron atoms and ions corresponds to the very large 
quantum numbers n of the hydrogen-like ions with the 
same Z. Therefore, the function F{aZ) in (§ cor¬ 
responds to the limit n, n' —> oo for the function 
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TABLE I: Calculated spectra of Mg I and A1 II. Experimental data are from Ref. [34]. Theoretical results are presented in 
columns 4-7: Cl is the calculation within configuration interaction method, MBPT - combined Cl and second order MBPT 
method, and AO - combined Cl and all-order method. Br+QED stands for contribution of the Breit interaction and QED. For 
the ground state we give binding energy of the two valence electrons (the sum of the first two ionization potentials). For other 
states we give excitation energies from the ground state. All energies are in cm“^. Next three columns give relative errors in 
percent for different methods. The Final theoretical values (CI-|-all-order-|-Breit-|-QED) and their errors are listed in the last 
three columns. For the final error we give both relative (in %) and absolute value (in cm“^). 


Element 

State 

Expt. 

Cl 

MBPT 

AO 

Br-|-QED Diff. with expt.(%) 
contrib. Cl MBPT AO 

Final theory 
AO+Br+QED 

Mg I 

3s^ 

'So 

182939 

179554 

182685 

182875 

-47 

1.9 

0.14 

0.035 

182828 

0.060 /111 


3s3p 

"Po 

21850 

20919 

21792 

21851 

-10 

4.3 

0.27 

-0.002 

21841 

0.042 / 9 


3s3p 


21871 

20939 

21814 

21872 

-11 

4.3 

0.26 

-0.009 

21862 

0.040 / 9 


3s3p 


21911 

20980 

21857 

21915 

-13 

4.2 

0.25 

-0.021 

21903 

0.037 / 8 


3s3p 

'Pi 

35051 

34471 

35030 

35052 

-15 

1.7 

0.06 

-0.002 

35037 

0.040 / 14 


3s4s 

"Si 

41197 

40419 

41132 

41184 

-14 

1.9 

0.16 

0.032 

41170 

0.066 / 27 


3s4s 

'So 

43503 

42678 

43438 

43490 

-14 

1.9 

0.15 

0.031 

43475 

0.064 / 28 


3s3d 

'P2 

46403 

45126 

46288 

46364 

-20 

2.8 

0.25 

0.085 

46344 

0.127 / 59 


3s4p 

"Po 

47841 

46932 

47766 

47822 

-16 

1.9 

0.16 

0.041 

47806 

0.073 / 35 


3s4p 

"Pi 

47844 

46935 

47769 

47825 

-16 

1.9 

0.16 

0.040 

47809 

0.073 / 35 


3s4p 

"P2 

47851 

46942 

47776 

47832 

-16 

1.9 

0.16 

0.039 

47816 

0.073 / 35 


3s3d 

"P2 

47957 

46972 

47851 

47904 

-18 

2.1 

0.22 

0.110 

47886 

0.147 / 71 


3s3d 

"Pa 

47957 

46972 

47851 

47904 

-18 

2.1 

0.22 

0.110 

47886 

0.147 / 71 


3s3d 

"Pi 

47957 

46972 

47851 

47904 

-18 

2.1 

0.22 

0.110 

47886 

0.147 / 71 


3s4p 

'Pi 

49347 

48497 

49284 

49332 

-16 

1.7 

0.13 

0.030 

49316 

0.063 / 31 

A1II 

3s2 

'So 

381308 

376617 

381166 

381313 

-103 

1.2 

0.04 

-0.001 

381210 

0.026 / 98 


3s3p 

"Po 

37393 

36257 

37354 

37406 

-14 

3.0 

0.10 

-0.035 

37392 

0.003 / 1 


3s3p 

"Pi 

37454 

36319 

37419 

37471 

-18 

3.0 

0.09 

-0.046 

37454 

0.000 / 0 


3s3p 

"P2 

37578 

36444 

37550 

37603 

-24 

3.0 

0.07 

-0.066 

37579 

-0.003 / -I 


3s3p 

'Pi 

59852 

59538 

59900 

59882 

-27 

0.5 

-0.08 

-0.050 

59855 

-0.005 /-3 


3p" 

'P2 

85481 

83516 

85415 

85494 

-45 

2.3 

0.08 

-0.015 

85450 

0.037 / 32 


3s4s 

"Si 

91275 

90010 

91238 

91288 

-32 

1.4 

0.04 

-0.015 

91256 

0.020 / 19 


3p" 

"Po 

94085 

92598 

94060 

94087 

-38 

1.6 

0.03 

-0.003 

94049 

0.038 / 36 


3p" 

"Pi 

94147 

92661 

94126 

94153 

-42 

1.6 

0.02 

-0.006 

94112 

0.038 / 36 


3p" 

"P2 

94269 

92783 

94255 

94283 

-48 

1.6 

0.01 

-0.015 

94234 

0.036 / 34 


3s4s 

'So 

95351 

94012 

95320 

95369 

-33 

1.4 

0.03 

-0.020 

95336 

0.015 / 15 


3s3d 

"Ps 

95549 

94171 

95443 

95463 

-45 

1.4 

0.11 

0.091 

95418 

0.138 /131 


3s3d 

"P2 

95551 

94172 

95444 

95464 

-45 

1.4 

0.11 

0.091 

95419 

0.138 /131 


3s3d 

"Pi 

95551 

94172 

95444 

95464 

-44 

1.4 

0.11 

0.091 

95420 

0.138 /132 


3s4p 

"Po 

105428 

103965 

105379 

105431 

-35 

1.4 

0.05 

-0.004 

105396 

0.030 / 32 


3s4p 

"Pi 

105442 

103979 

105393 

105446 

-36 

1.4 

0.05 

-0.004 

105410 

0.030 / 31 


3s4p 

"P2 

105471 

104009 

105424 

105477 

-38 

1.4 

0.04 

-0.006 

105440 

0.030 / 31 


3s4p 

'Pi 

106921 

105578 

106885 

106929 

-37 

1.3 

0.03 

-0.008 

106892 

0.027 / 28 


Fn^n'{c(Z). According to [30] this function slowly de¬ 
creases with Z, but for the range 12 < Z < 17 we can 
take F{aZ) « 4. Thus, we can write: 

^ o./o „3 v2 fc\ 

~ 8v2 a Z . (6) 

7r(Zi -I- Ij 

In our CI-I-AO calculations we add these QED corrections 
to the one-electron radial integrals of valence electrons. 

III. RESULTS AND DISCUSSIONS 

In Tables 0-1111 we compare with the experiment the 
results of our Cl, CI-I-MBPT, and CI-I-AO ab initio calcu¬ 
lations of the magnesium isoelectronic sequence. Table 0 


includes ions Mg I, A1 II; Table jTT] and Table [Till present 
results for Si III and P IV and for S V and Cl VI respec¬ 
tively. For the ground states we give the two-electron 
binding energies. These energies are equal to the sum of 
the first two ionization potentials. For all other states 
we give transition energies from the ground states. All 
energy values are in cm“^. The same notations are used 
in all tables. To illustrate the accuracy of each of the the¬ 
oretical approaches the relative differences of our results 
with the experiment are given in the columns designated 
as “Diff. with expt.”. 

Our main concern here is the analysis of the accuracy 
of the three theoretical methods and the role of different 
corrections. On Fig.j^and Fig.[^we present relative the¬ 
oretical errors for different methods for the four typical 
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TABLE II: Same as in Table U for Si III and P IV . 


Element 

State 

Expt. 

Cl 

MBPT 

AO 

Br-tQED 

Diff. 

with expt.(%) 

Final theory 









contrib. 

Cl 

MBPT 

AO 

AO+Br+QED 


Si III 

3s^ 

^So 

634232 

628587 

634174 

634298 

-181 

0.9 

0.01 

-0.010 

634117 

0.018 

/ 

115 


3s3p 

^Po 

52725 

51546 

52702 

52750 

-18 

2.2 

0.04 

-0.049 

52733 

-0.015 

/ 

-8 


3s3p 

^Pi 

52853 

51678 

52838 

52887 

-24 

2.2 

0.03 

-0.063 

52863 

-0.018 

/ 

-9 


3s3p 

^P2 

53115 

51944 

53114 

53163 

-37 

2.2 

0.00 

-0.090 

53126 

-0.021 

/ 

-11 


3s3p 

^Pl 

82884 

83001 

82978 

82942 

-39 

-0.1 

-0.11 

-0.069 

82903 

-0.023 

/ 

-19 


3p" 

^D2 

122215 

120238 

122216 

122286 

-69 

1.6 

0.00 

-0.059 

122217 

-0.002 

/ 

-2 


3p" 

"Po 

129708 

128580 

129724 

129731 

-52 

0.9 

-0.01 

-0.018 

129679 

0.022 

/ 

29 


3p" 


129842 

128715 

129864 

129872 

-59 

0.9 

-0.02 

-0.023 

129813 

0.022 

/ 

29 


3p" 

®P2 

130101 

128977 

130138 

130146 

-72 

0.9 

-0.03 

-0.035 

130075 

0.020 

/ 

26 


3s3d 


142944 

141714 

142867 

142854 

-81 

0.9 

0.05 

0.063 

142773 

0.119 

/ 

170 


3s3d 


142946 

141715 

142867 

142854 

-79 

0.9 

0.05 

0.064 

142775 

0.119 

/ 

171 


3s3d 


142948 

141716 

142868 

142855 

-78 

0.9 

0.06 

0.065 

142777 

0.120 

/ 

172 


3s4s 

"Si 

153377 

151780 

153378 

153425 

-57 

1.0 

0.00 

-0.031 

153368 

0.006 

/ 

9 


3p" 

^So 

153444 

152676 

153628 

153610 

-70 

0.5 

-0.12 

-0.108 

153540 

-0.062 

/ 

-96 


3s4s 

^So 

159070 

157563 

159096 

159136 

-60 

0.9 

-0.02 

-0.042 

159076 

-0.004 

/ 

-7 


3s3d 

^D2 

165765 

165093 

165825 

165778 

-75 

0.4 

-0.04 

-0.008 

165703 

0.037 

/ 

62 


3s4p 

"Po 

175230 

173427 

175215 

175264 

-61 

1.0 

0.01 

-0.019 

175203 

0.016 

/ 

27 


3s4p 

"Pi 

175263 

173461 

175250 

175299 

-63 

1.0 

0.01 

-0.021 

175236 

0.015 

/ 

27 


3s4p 

"P2 

175336 

173536 

175327 

175376 

-66 

1.0 

0.01 

-0.023 

175310 

0.015 

/ 

26 


3s4p 

"Pi 

176487 

174829 

176487 

176530 

-65 

0.9 

0.00 

-0.024 

176465 

0.013 

/ 

22 

P IV 

3s^ 

"So 

939386 

933026 

939426 

939530 

-283 

0.7 

0.00 

-0.015 

939247 

0.015 

/ 

139 


3s3p 

"Po 

67918 

66765 

67907 

67951 

-20 

1.7 

0.02 

-0.048 

67931 

-0.019 

/ 

-13 


3s3p 

"Pi 

68146 

67001 

68149 

68192 

-30 

1.7 

0.00 

-0.068 

68162 

-0.023 

/ 

-16 


3s3p 

"P2 

68615 

67479 

68641 

68686 

-52 

1.7 

-0.04 

-0.103 

68634 

-0.027 

/ 

-19 


3s3p 

"Pi 

105190 

105754 

105320 

105277 

-52 

-0.5 

-0.12 

-0.082 

105225 

-0.033 

/ 

-34 


3p^ 

"P2 

158142 

156361 

158207 

158266 

-96 

1.1 

-0.04 

-0.079 

158170 

-0.018 

/ 

-28 


3p2 

"Po 

164941 

164237 

164987 

164984 

-66 

0.4 

-0.03 

-0.026 

164918 

0.014 

/ 

23 


3p^ 

"Pi 

165185 

164484 

165243 

165240 

-77 

0.4 

-0.04 

-0.033 

165163 

0.014 

/ 

23 


3p" 

"P2 

165654 

164962 

165739 

165736 

-100 

0.4 

-0.05 

-0.049 

165636 

0.011 

/ 

18 


3s3d 

"P2 

189398 

188560 

189368 

189335 

-125 

0.4 

0.02 

0.033 

189210 

0.099 

/ 

188 


3s3d 

"Pi 

189398 

188557 

189363 

189330 

-121 

0.4 

0.02 

0.036 

189209 

0.100 

/ 

189 


3s3d 

"Ps 

189398 

188555 

189360 

189327 

-119 

0.4 

0.02 

0.037 

189209 

0.100 

/ 

189 


3p2 

"So 

194592 

194530 

194857 

194821 

-96 

0.0 

-0.14 

-0.118 

194725 

-0.068 

/- 

-133 


3s3d 

"P2 

219154 

219299 

219300 

219215 

-116 

-0.1 

-0.07 

-0.028 

219099 

0.025 

/ 

55 


3s4s 

"Si 

226898 

225066 

226947 

226990 

-89 

0.8 

-0.02 

-0.040 

226901 

-0.001 

/ 

-3 


3s4s 

"So 

233998 

232162 

234071 

234110 

-92 

0.8 

-0.03 

-0.048 

234018 

-0.009 

/ 

-20 


3s4p 

"Po 

256553 

254539 

256585 

256629 

-94 

0.8 

-0.01 

-0.029 

256535 

0.007 

/ 

19 


3s4p 

"Pi 

256612 

254602 

256647 

256691 

-96 

0.8 

-0.01 

-0.031 

256594 

0.007 

/ 

18 


3s4p 

"P2 

256760 

254753 

256803 

256847 

-104 

0.8 

-0.02 

-0.034 

256743 

0.007 

/ 

17 


3s4p 

"Pi 

257523 

255651 

257570 

257610 

-102 

0.7 

-0.02 

-0.034 

257508 

0.006 

/ 

14 


ions of the sequence including lightest and heaviest ones. 
For the ground states we again give the errors for the 
two-electron binding energies. For other states the errors 
correspond to the transition frequencies from the ground 
states. The Plots in Fig. [T] demonstrate the accuracy of 
all three methods for the Dirac-Coulomb Hamiltonian. 
We see that for each method there is certain improve¬ 
ment along the isoelectronic sequence. For example, the 
average accuracy of the Cl method for Mg I is about 2%, 
and improves to roughly 0.4% for the Cl VI. Similarly the 
accuracy of the CI-I-MBPT method improves from 0.2% 
for Mg I to 0.1% for Cl VI. At the same time the dif¬ 
ference between CI-I-MBPT and CI-fAO decreases with 
the ion charge Z, and almost disappears for Cl VI. This 


indicates smaller role of the higher order core-valence cor¬ 
relations for heavier ions. It is interested to note that the 
Cl space for the triplet states was already saturated on 
the 15spdfg level. But in order to obtain similar accu¬ 
racy for singlet states we had to increase the Cl space to 
2lspdfg. 

Fig. [^demonstrates the influence of the Breit and QED 
corrections. These corrections are absolutely negligible 
for the Cl calculations and only marginally noticeable 
for the CH-MBPT method. Therefore, Fig. presents 
results only for the CH-AO method. We see that for the 
light ions (Mg I, A1 II) these corrections are negligible 
at the existing level of accuracy of treating electron cor¬ 
relations. However, for the heaviest ion (Cl VI) these 
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TABLE III: Same as in Table |I] for S V and Cl VI . 


Element 

State 

Expt. 

Cl 

MBPT 

AO 

Br-bQED 

contrib. 

Diff. 

Cl 

with expt.(%) 
MBPT AO 

Final theory 
AO+Br+QED 


S V 


'5o 

1295709 

1288803 

1295858 

1295948 

-415 

0.5 

-0.01 

-0.018 

1295532 

0.014 

/ 

176 


3s3p 


83024 

81926 

83020 

83060 

-21 

1.3 

0.00 

-0.044 

83039 

-0.018 

/ 

-15 


3s3p 

®Pi 

83394 

82307 

83409 

83449 

-37 

1.3 

-0.02 

-0.067 

83412 

-0.022 

/ 

-18 


3s3p 

®P2 

84155 

83087 

84208 

84248 

-70 

1.3 

-0.06 

-0.111 

84178 

-0.027 

/ 

-23 


3s3p 

^Pl 

127151 

128134 

127312 

127265 

-67 

-0.8 

-0.13 

-0.090 

127198 

-0.037 

/ 

-48 


3p" 

^D2 

193739 

192223 

193859 

193910 

-126 

0.8 

-0.06 

-0.088 

193784 

-0.023 

/ 

-45 


3p" 

"Po 

199967 

199690 

200036 

200027 

-80 

0.1 

-0.03 

-0.030 

199947 

0.010 

/ 

20 


3p" 

"Pi 

200371 

200101 

200460 

200451 

-99 

0.1 

-0.04 

-0.040 

200352 

0.009 

/ 

18 


3p" 

"P2 

201146 

200891 

201277 

201268 

-134 

0.1 

-0.07 

-0.061 

201135 

0.006 

/ 

11 


3s3d 

"Pi 

234942 

234556 

234947 

234905 

-166 

0.2 

0.00 

0.016 

234739 

0.086 

/ 

203 


3s3d 

"P2 

234947 

234566 

234959 

234917 

-172 

0.2 

-0.01 

0.013 

234746 

0.086 

/ 

201 


3s3d 

"Ps 

234956 

234580 

234977 

234936 

-180 

0.2 

-0.01 

0.009 

234756 

0.085 

/ 

200 


3p" 


235350 

235873 

235675 

235633 

-126 

-0.2 

-0.14 

-0.120 

235507 

-0.067 

/- 

-157 


3s3d 

"P2 

270700 

271683 

270933 

270829 

-165 

-0.4 

-0.09 

-0.048 

270664 

0.013 

/ 

36 


3s4s 

"5i 

311595 

309588 

311692 

311730 

-129 

0.6 

-0.03 

-0.043 

311601 

-0.002 

/ 

-6 


3s4s 

"So 

320108 

318110 

320238 

320274 

-134 

0.6 

-0.04 

-0.052 

320140 

-0.010 

/ 

-32 


3p3d 

"P2 

323133 

321355 

323119 

323130 

-200 

0.6 

0.00 

0.001 

322931 

0.062 

/ 

202 


3p3d 

"Ps 

323547 

321780 

323559 

323571 

-224 

0.5 

0.00 

-0.007 

323348 

0.062 

/ 

200 


3p3d 

"P4 

324080 

322331 

324126 

324139 

-254 

0.5 

-0.01 

-0.018 

323885 

0.060 

/ 

195 


3p3d 

"P2 

328454 

327374 

328614 

328591 

-234 

0.3 

-0.05 

-0.042 

328357 

0.030 

/ 

97 


3p3d 

"P2 

345338 

344679 

345437 

345389 

-225 

0.2 

-0.03 

-0.015 

345165 

0.050 

/ 

173 


3p3d 

"Pi 

345713 

345038 

345818 

345771 

-232 

0.2 

-0.03 

-0.017 

345539 

0.050 

/ 

174 


3p3d 

"Po 

345953 

345255 

346062 

346016 

-241 

0.2 

-0.03 

-0.018 

345775 

0.051 

/ 

178 

Cl IV 

3s^ 

"So 

1702996 

1695325 

1702922 

1703001 

-580 

0.5 

0.00 

0.000 

1702421 

0.034 

/ 

575 


3s3p 

"Po 

98062 

97047 

98079 

98115 

-22 

1.0 

-0.02 

-0.054 

98093 

-0.032 

/ 

-31 


3s3p 

"Pi 

98621 

97623 

98664 

98701 

-44 

1.0 

-0.04 

-0.082 

98657 

-0.036 

/ 

-36 


3s3p 

"P2 

99782 

98813 

99879 

99916 

-92 

1.0 

-0.10 

-0.134 

99824 

-0.042 

/ 

-42 


3s3p 

"Pi 

148947 

150326 

149149 

149102 

-85 

-0.9 

-0.14 

-0.104 

149018 

-0.047 

/ 

-71 


3p2 

"P2 

229219 

228003 

229404 

229448 

-160 

0.5 

-0.08 

-0.100 

229288 

-0.030 

/ 

-69 


3p" 

"Po 

234886 

235038 

234996 

234983 

-94 

-0.1 

-0.05 

-0.041 

234888 

-0.001 

/ 

-2 


3p2 

"Pi 

235518 

235681 

235655 

235643 

-123 

-0.1 

-0.06 

-0.053 

235520 

-0.001 

/ 

-2 


3p" 

"P2 

236721 

236907 

236922 

236911 

-175 

-0.1 

-0.08 

-0.080 

236735 

-0.006 

/ 

-14 


3p2 

"So 

275988 

277074 

276389 

276345 

-162 

-0.4 

-0.15 

-0.129 

276184 

-0.071 

/- 

-196 


3s3d 

"Pi 

279758 

279851 

279830 

279787 

-221 

0.0 

-0.03 

-0.010 

279566 

0.069 

/ 

192 


3s3d 

"P2 

279773 

279877 

279860 

279817 

-230 

0.0 

-0.03 

-0.016 

279587 

0.067 

/ 

186 


3s3d 

"Ps 

279804 

279916 

279906 

279863 

-245 

0.0 

-0.04 

-0.021 

279618 

0.066 

/ 

186 


3s3d 

"P2 

320925 

322706 

321257 

321147 

-222 

-0.6 

-0.10 

-0.069 

320925 

0.000 

/ 

0 


3p3d 

"P2 

383627 

382362 

383715 

383722 

-258 

0.3 

-0.02 

-0.025 

383464 

0.042 

/ 

163 


3p3d 

"Ps 

384280 

383037 

384412 

384420 

-296 

0.3 

-0.03 

-0.036 

384124 

0.041 

/ 

156 


3p3d 

"P4 

385115 

383896 

385294 

385302 

-342 

0.3 

-0.05 

-0.049 

384960 

0.040 

/ 

155 


3p3d 

"P2 

389464 

388983 

389717 

389687 

-312 

0.1 

-0.06 

-0.057 

389376 

0.023 

/ 

88 


3s4s 

"Si 

407323 

405211 

407496 

407531 

-180 

0.5 

-0.04 

-0.051 

407351 

-0.007 

/ 

-28 


3p3d 

"P2 

409004 

409106 

409196 

409142 

-297 

0.0 

-0.05 

-0.034 

408845 

0.039 

/ 

159 


3p3d 

"Pi 

409551 

409657 

409750 

409696 

-304 

0.0 

-0.05 

-0.035 

409392 

0.039 

/ 

159 


3p3d 

"Po 

409949 

410028 

410153 

410101 

-319 

0.0 

-0.05 

-0.037 

409783 

0.041 

/ 

166 


corrections become essential and somewhat improve the 
final accuracy: an average error decreases from 0.06% to 
0.04%. Note that for the ions considered here Breit and 
QED corrections are comparable. Breit interaction gen¬ 
erally improves the fine structure splittings, while QED 
corrections decrease the overall scatter of errors. 

Einally, we consider the energy dependence of the ef¬ 
fective Hamiltonian. We find out that respective correc¬ 
tions are comparable in size to the difference between the 
CI-I-MBPT and CI-I-AO methods. It agrees with the con¬ 


clusion in Ref. [11] that accurate treatment of the high- 
order corrections requires also including corrections on 
the energy dependence. The average size of these correc¬ 
tions to the valence energies monotonously decrease from 
0.06% for Mg I to 0.02% for Cl VI. In general they do 
not improve the agreement with the experiment. Only for 
the CI-I-MBPT calculations of Mg I some improvement 
(about 0.1%) take place. This may mean that correc¬ 
tions on the energy dependence cancel some high-order 
terms which are missing in our calculations. We con- 
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FIG. 1: The relative theoretical errors for Mg-like ions within Cl, CI+MBPT, and CI+AO methods for the Dirac-Conlomb 
Hamiltonian. 


elude that for the present variant of the CI+AO method, 
which is based on the linearized SD CC, these corrections 
should be neglected. They also should not be included in 
the CI+MBPT calculations. Note that this significantly 
simplihes calculations with the package [32]. 


IV. CONCLUSION 

In this paper we studied the accuracy of the CI+AO 
method la HD for the isoelectronic sequence of Mg. 
These ions have ten electrons in the closed shells and two 
valence electrons and are often used as a test ground for 
the atomic theory. We found out that CI+AO method 
provides higher accuracy than the simpler and more com¬ 
mon CI+MBPT method [1]. While the accuracy of the 
CI+MBPT method was on the level 0.1 - 0.2%, the accu¬ 
racy of the CI+AO was roughly two times higher, 0.05 - 
0.1%. Note that conventional valence Cl is an order of 
magnitude less accurate. 

The accuracy slightly increases along the isoelectronic 
sequence. For the first member of the sequence, Mg I, 
the final accuracy of the theory for the low-lying levels 
is close to 0.08% and for the last ion. Cl VI, it is about 


0.04%. Breit and QED corrections start to become im¬ 
portant on this level of accuracy for the atoms and ions 
with Z > 20. Retardation part of the Breit interaction 
is known to be signihcantly smaller than magnetic part 
and can be still neglected. For QED corrections it is suf¬ 
ficient to account only for the s-wave contribution and 
use simplified semiempirical expression Q. 

We also studied corrections on the energy dependence 
of the effective Hamiltonian in the CI+AO method. On 
the one hand, we found them to be rather small. On the 
other hand, these corrections did not improve agreement 
with the experiment. We conclude that corrections on 
the energy dependence can be neglected for the present 
variant of the CI+AO method, when the all-order part 
corresponds to the linearized coupled cluster method in 
the SD approximation. This significantly simplihes cal¬ 
culations and makes the whole method more practical. 
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